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For describing the first-order isostructural valence phase 
transition in mixed valence compounds we develop a new ap- 
proach based on the lattice Anderson model. We take into 
account the Coulomb interaction between localized / and 
conduction band electrons and two mechanisms of electron- 
lattice coupling. One is related to the volume dependence of 
the hybridization. The other is related to local deformations 
produced by / shell size fluctuations accompanying valence 
fluctuations. The large /-state degeneracy allows us to use 
the 1/N expansion method. Within the model we develop 
a mean-field theory for the first-order valence phase transi- 
tion in YbInCu4. ft is shown that the Coulomb interaction 
enhances the exchange interaction between / and conduction 
band electron spins and is the driving force of the phase tran- 
sition. A comparison between the theoretical calculations and 
experimental measurements of the valence change, suscep- 
tibility, specific heat, entropy, elastic constants and volume 
change in YbInCu4 and YbAgCu4 are presented, and a good 
quantitative agreement is found. On the basis of the model 
we describe the evolution from the first-order valence phase 
transition to the continuous transition into the heavy-fermion 
ground state in the series of compounds YbIni_ a; Ag I Cu4. 
The effect of pressure on physical properties of YblnCu4 is 
studied and the H — T phase diagram is found. 
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I. INTRODUCTION 

The first-order isostuctural valence phase transition 
observed in rare earth intermetallic compounds is one of 
the most striking and interesting phenomena in strongly 
correlated electron systems. For the first time the 
transition was discovered in elemental Ce which at ap- 
plied pressure shows a a — > 7 transition (for a review, 
see Other compounds demonstrating the transi- 

tion are samarium monochalcogenides, |a,0j a few Eu- 
based compounds, |||(| CeNi^Coi-^Sn, [£ andYbInCu4. 
|s| |Tof The latter compound has attracted much recent 
attention because YbInCu4 demonstrates a first-order 
isostuctural valence phase transition at ambient pressure. 
The remarkable peculiarity of the transition lies in the 
small change of valence of the Yb ions from nearly +3 at 
high temperatures (T >50 K) to about 2.8-2.9 below the 



critical temperature T v — 42 K. The volume expansion 
of YblnCu4 taking place at T = T v is also very small, 
about +0.5%, whereas in Ce the volume change achieves 
-15%. Despite the small valence and volume changes 
at T = T v , many physical properties of YbInCu4 are 
changed dramatically. Already first magnetic measure- 
ments J9|,|l0| showed that at the transition the susceptibil- 
ity drops abruptly, and the Curie- Weiss temperature de- 
pendence above T v is replaced by an almost temperature 
independent behavior which is typical for a nonmagnetic 
metal with an enhanced Pauli paramagnetic susceptibil- 
ity. While the Curie- Weiss behavior at high temperatures 
can be easily related to the magnetic moment J = 7/2 
of trivalent Yb ions, the formation of the Pauli paramag- 
netic state is less clear, taking into account that the va- 
lence change ANf at T v is estimated to be only about 0.1 
as it was revealed by x-ray absorption measurements at 
the Yb L3 edge. f|,[llj From the same point of view the re- 
sults of Hall effect measurements |ll],[l2| look very amaz- 
ing, since they give evidence for a very large change of the 
charge carrier concentration from 0.07 per formula unit 
at high temperatures to 2.2 per formula unit below T v . In 
the high temperature phase the Hall measurements are 
consistent with the band calculations |ll| which showed 
that in this phase YbInCu4 has a semimetal band struc- 
ture similar to LuInCu4. Transport measurements in 
YbInCu4 also shows that the compound is semimetalic 
above T v . [|l4| Resistance measurements revealed a large 
and abrupt drop of electrical resistivity when decreas- 
ing the temperature below T v . Investigations of Cu nu- 
clear quadrupolar resonance |]l5|-p7| (NQR) and Cu and 
In Knight shift in YbInCu 4 show that at T > T v 
the 4/ electrons are localized and interact weakly with 
the conduction band electrons while at T < T v a Fermi- 
liquid state with delocalized Yb 4/ electrons is formed. 
The measurements revealed occurrence of a strong hy- 
bridization between Yb 4/ electron states and conduc- 
tion band states in the Fermi-liquid state. This result is 
in agreement with Mossbauer experiments |J that also 
revealed a nonmagnetic state of the Yb ions below T v . 
Additional evidence for the fact comes from the temper- 
ature behavior of the entropy (S) inferred from specific 
heat measurements. Above T v the entropy is large 
and can be explained by a large contribution from the 
/ electrons. At T v the entropy drops significantly and 
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with decreasing temperature it tends to zero in the limit 
T — > 0. It means that at T = the ground state of 
the electron system consisting of conduction band and 
/ electrons, is a singlet state which in accordance with 
magnetic, NQR and transport measurements is a Fermi- 
liquid state. 

Summarizing, one can conclude that the main mystery 
of YbInCu4 is in the following. In spite of a small valence 
change at the phase transition (ANf ~ 0.1) from which 
one could expect that / electrons would save their local- 
ized character, the thermal, magnetic, NQR and trans- 
port measurements mentioned above clearly demonstrate 
that below T v the / electrons loose their localized charac- 
ter and form a nonmagnetic Fermi-liquid state together 
with the conduction electrons . 

Taking into account the nonmagnetic character of the 
Yb ions at T < T v , the enhanced values of the zero 
temperature susceptibility y(0) and linear specific heat 
coefficient 7, Felner et al. || came to the conclusion that 
YbInCu4 is a "light heavy-fermion system" . 

The uniqueness of the first-order valence phase transi- 
tion in YbInCu4 becomes especially clear when we com- 
pare the compound with other isostructural Yb com- 
pounds, namely, YbXCu 4 with X=Ag, Au, Cd, Mg, Tl, 
and Zn. A detailed study and comparison of the com- 
pounds have recently been performed by Sarrao et al. 
p0[ Though the compounds have the same face-centered- 
cubic crystal structure, their physical properties are very 
different. Among the compounds YbAgCu 4 is of spe- 
cial interest. It is a typical heavy-fermion (HF) com- 
pound with no magnetic order. A substitution of Ag 
by In results in crossover from the HF behavior to the 
first-order isostructural valence phase transition, as it has 
been found when investigating the series of compounds 
Ybln^Ag^C^. 

In order to explain the valence phase transition in Ce 
several models have been proposed. In the framework of 
promotional models [23|j24 | it is assumed that there is an 
explicit change of the electronic configuration of Ce ions 
from 4/ 1 to 4/°, and / electrons transfer into the con- 
duction band. A Mott transition picture considers the 
valence transition as a localization-delocalization transi- 



tion in the system of / electrons. 25 The Kondo volume- 
collapse (KVC) model p^^7| explains the valence phase 
transition by a strong volume dependence of the Kondo 
temperature 2\, neglecting correlation effects between 
the Kondo effect on different / ions. Each of the mod- 
els has advantages and shortcomings. The shortcomings 
becomes especially contrasting when we try to use the 
models for explaining an evolution of the valence phase 
transition in the series of compounds YbIni_ :c Ag 2: Cu4. 

The promotional models demand the valence change 
to be about 1, while in YbIni_ :c Ag :c Cu4 the change is 
only about 0.1. This disagreement does not allow within 
the models to explain quantitatively, using realistic phys- 
ical parameters, magnetic and other properties of the Yb 
compounds in the low temperature phase. A recent at- 
tempt 1 28 1 to correct the deficiency of the Falicov-Kimbal 



model is based on the assumption that only a fraction 
of Yb ions are active and other ones are strongly hy- 
bridized. However, the assumption has no experimental 
confirmation. It contradicts , for example, recent ther- 
moanalytical investigations [^9| of the system Yb-In-Cu 
which showed that by choosing certain starting composi- 
tions one can obtain such YbInCu4-samples which have 
T v ~ 40 K, an improved site order and a reduced proba- 
bility of Yb ions to occupy In-sites. The Falicov-Kimbal 
model assumes that the 4/ 14 configuration is the ground 
state of Yb ions and 4/ 13 configuration is the first excited 
state. This assumption contradicts with a large number 
of experimental data that give direct and indirect evi- 
dences for the fact that the 4/ 13 state is the ground state, 
and the 4/ 13 configuration is the first excited state (see 
below Sec. II). The promotional model actually neglects 
a hybridization between / states and conduction band 
states. It contradicts with experimental data on NQR 
]l5| [nj and the Knight shift measurements 18 men- 
tioned above. Moreover, on the basis of the models it 
is impossible to explain the development of a continuous 
transition into the HF ground state in Yblni^AgzCv^ 
at x > 0.2 p2fl , as the Kondo effect taking place in the 
compounds, can not be described in the framework of the 
model. 

In terms of the Mott transition picture [^| at T < T v 
the / band width W has to become larger than the on- 
site Coulomb repulsion energy U between / electrons. It 
is possible if the volume of the lattice shrinks below T v . 
Such a situation takes place in Ce, but not in YbInCu 4 
. The volume of YbInCu4 increases by 0.5%, which evi- 
dences against the Mott transition picture. 

It should be mentioned that both the promotional 
models and Mott transition picture have a charge fluc- 
tuation energy as the energy scale. Usually a charge 
fluctuation energy is too large to provide a convinc- 
ing explanation of the small energy scales observed in 
YbIni_ :c Ag 2; Cu4. 

The KVC model faces a serious problem when explain- 
ing pressure dependence of the zero temperature suscep- 
tibility in YbInCu4. The model needs the Gruneisen pa- 
rameter r to be negative and equal to —3000, an un- 
physical value, while the experiment jl9| gives T rs —31. 
The negative sign of T also constitutes a challenge to the 
model because within the usual mechanism of the volume 
dependence of the Kondo temperature due to the vol- 
ume dependence of the hybridization, T should be pos- 
itive. Moreover, the KVC model neglects correlations, 
or in other words, coherence in the Kondo screening of 
nearby Yb ions in the low temperature phase. However, 
it is the latter effect that results in the formation of the 
heavy fermion state in a heavy fermion compound such 
as YbAgCu4(see, for example, reviews [ j30"| , fn"| ) . 

In the present paper we propose a new theoretical ap- 
proach that allows us both to obtain quantitative agre- 
ment with various experimental data and to avoid con- 
tradictions and deficiencies of the previous approaches. 
Our approach is based on the conclusion obtained from 
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the experimental data mentioned above that the low tem- 
perature phase of YbIni_ 2 .Ag a; Cu4 is the HF state. From 
this point of view the first-order valence phase transition 
is a normal metal (semimetal)-HF metal transition. Be- 
sides, it is supposed that at T > T v in the normal state 
the / electrons are localized on the / shells and produce 
the magnetic moments J = 7/2 of the Yb ions. In the 
HF state, i.e. at T < T v , f states are strongly hybridized 
with conduction band states, and as a result, quasiparti- 
cles near the Fermi surface have a hybrid character and 
an enhanced mass. Such a physical picture can be suc- 
cessfully described on the basis of the lattice Anderson 
model. However, it is important to note that the con- 
ventional lattice Anderson model allows to study only a 
continuous transition from the normal state with the in- 
coherent Kondo scattering into the heavy-fermion state 
with the coherent Kondo effect, as it occurs in heavy- 
fermion compounds. [ jSOfl In the present paper we extend 
the model in the following way. At first, we take into 
account a Coulomb repulsion between / and conduction 
band electrons. Second, we take into account two mech- 
anisms of the electron-lattice coupling. The first mech- 
anism is related to the volume dependence of the hy- 
bridization between / and conduction band states due to 
a change of the overlapping between wave functions of the 
conduction band and / electron states. It is the conven- 
tional mechanism of the electron-lattice interaction in HF 
compounds, that has been used specifically in the KVC 
model. However, it is necessary to note that valence fluc- 
tuations of the Yb ions are accompanied by fluctuations 
of the / shell size, which results in local lattice defor- 
mations. It is the second (/-shell-size-fluctuation) mech- 
anism of the electron lattice coupling considered in our 
approach. If one neglects crystal field splitting, then the 
/ level has the degeneracy N = 2 J + 1 = 8. Such a large 
degeneracy enables us to use the 1/N expansion tech- 
nique for attacking the extended lattice Anderson model. 
In the present work we shall calculate the free energy in 
the mean-field approximation, that corresponds to the 
leading order in 1/N. A calculation of 1/N corrections 
is beyond the scope of our paper. 

Our paper has the following structure. In Sec. II we in- 
troduce the extended lattice Anderson model. The mean- 
field approach for studying the model is developed in Sec. 
III. The electron-lattice coupling and volume change are 
considered in Sec. IV. In the next Sec. V we investigate 
magnetic properties of the system and magnetic field ef- 
fect. Results of numerical calculations based on the mean 
field approach are discussed in Sec. VI where we com- 
pare the theoretical results on temperature and pressure 
dependences of thermal, magnetic and elastic properties 
of YbInCu4 and YbAgCu4 with experimental measure- 
ments. In Sec. VII we will discuss the role of fluctuations 
around the mean-field solution in low and high tempera- 
ture phases and make conclusions.. 



II. MODEL 

A number of experimental data give direct evidences 
for the fact that Yb ions in YbInCu4 fluctuate between 
the magnetic ground state Yb +3 having the configura- 
tion 4/ 13 with one hole in the / shell and total magnetic 
moment j=7/2, and the excited nonmagnetic state Yb +2 
having a completely filled / shell with the configuration 
4/ 14 . One can mention the Mossbauer spectroscop y || 
and the photoemission spectra (PES) measurements | p2| . 
Moreover, according to the PES measurements, neither 
4/ 12 nor 4/ 13 2 F 5 / 2 contribute to low energy properties. 
Therefore, one can neglect transitions into these states. 

For describing the system it is convenient to use a hole 
representation. Valence fluctuations of Yb ions can be 
described in terms of hole transitions between / shells 
and the hole conduction band. The picture of behavior 
of an 4/ ion in a metal was developed by Hirst |3j|. It 
underlies the theoretical model we will use in the paper. 
The system of interacting / holes and conduction band 
holes can be described using the lattice Anderson Hamil- 
tonian in the slave boson representation: |30| , p4] , ^5|| 

Ha=J2 £ k c ^k c "k + E f f£n fan 
ak an 

+ J2( Vb t4nf a n + H.C.), (1) 

an 

where c^ k and c Q k are creation and annihilation opera- 
tors of conduction band holes with z-component a of the 
total magnetic moment, a = + l,---j, and wave 

vector k. Operators f£ n and f an create and annihilate 
/ holes on a lattice site with a number n. The Bose op- 
erators b„ and b n are the conventional slave bosons. A 
wave function of the state with one boson on a site n is 
equal to 6+ | O > . This state corresponds to the 4/ 14 
configuration of the Yb +2 ion on the site n. The operator 
b+b n determines the probability to find the / ion in the 
divalent state. In order to admit valent fluctuations only 
between Yb +3 and Yb +2 states the following constraint 
is imposed on each lattice site: 

b+b n +J2fanfan = l. (2) 

a 

We should mention that in the limit \i — Ef << ir \ V\ 2 po 
(p is the chemical potential and p$ is a conduction band 
density of states) one can integrate over the bosons. 
As a result, this model becomes the SU(N) periodic 
Coqblin-Schriffer model with an exchange interaction 
J =\V\ 2 /{p-E f ): @@ 

H A=Y1 £ k4k C «k ~ J ^2(Ca"nfan)(fp n Cf3 n ) (3) 
ak af3n 

This model describes the integral valence case when the 
f 13 state behaves like a spin interacting antiferromag- 
netically with conduction band holes. In mixed-valence 
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compounds Ef can be close to [i, i.e. /i — Ef ~ 7r |U| po- 
ifjlf Therefore, in the latter case it is necessary to use the 
lattice Anderson model. 

When a / hole transits into the conduction band, 
then a corresponding Yb ion acquires an excess negative 
charge (Sq = e). The excess charge produces a Coulomb 
potential that attracts conduction band holes. Due to the 
small charge carrier concentration (Nq c = 0.07 per for- 
mula unit) the interaction is not completely screened. At 
least, the Coulomb interaction between the fluctuating 
charge of Yb ions and electrons (or holes) on neighboring 
Cu ions of which p and d states mainly contribute to the 
formation of the conduction band, should be taken into 
account. Unfortunately, a detailed consideration of the 
long range Coulomb interaction is a complicated prob- 
lem. For the sake of simplicity we approximate the in- 
teraction by use of an on-site Coulomb attraction be- 
tween the excess charge and conduction holes. The same 
approximation is made within the promotional models. 
p3| , p4j In the slave boson representation this on-site at- 
traction can be written as follows: 



U c f ^ ~] bu^n^a ,, ' ■mi - 



(4) 



where the parameter U c f > 0. It characterizes the ef- 
fective value of the Coulomb interaction between 4/ and 
conduction band holes. Using the constraint Eq.(j|), one 
can rewrite the interaction in the form 



a(3n 



u, 



cf 



an 



(5) 



The first term describes an on-site Coulomb repulsion 
between conduction band and / holes. The second term 
results in an energy shift of the conduction band. Accord- 
ing to the representation (^), the attractive interaction 
Eq. (Q) is equivalent to the on-site repulsion between con- 
duction band and / electrons. In the present paper we 
prefer to use the representation Eq.(Q). 

The total Hamiltonian of the system under considera- 
tion is given by 



H = H A + H C 



(6) 



The model described by the Hamiltonian can be solved 
exactly in the limit of large degeneracy N = 2j + 1 >> 1 
if we suppose the following dependence of the parameters 
V and U cf on N: 



V = VN~ 1/2 , Ucf = UcfN" 



(7) 



where V and U c f are finite in the limit N » 1. More- 
over, it has to be supposed that the initial number of 
conduction band holes, tiq c = Nq c /N , and the number 
of / holes, rif)f — N f/N, per unit cell and orbital are 
finite in the limit N » 1. The constrain Eq.(^) must 
be replaced by the following one: |p6[ 



btbr. 



y ] fan fan 



Nn, 



of 



(8) 



where no/ = 1/8 for the considered Yb system. 

The partition function of the model Eq.(^|) may be 
written as a path integral over Grassmann variables c + , 
c, / + ,/ and Bose variables b + , b and A : 





Z = J D(c + cf + fb+b\)eM- J cLtL{t)). (9) 

o 

Here the Lagrangian L is given by 

L{t) = ^2 c tk( 9 r + £k - ")c«k + ^2 b+d T b n 

ak n 

+ Y.^n{9r + Ef - U)/ Q „ + {Vbtc+J an 

cm 

+H.c.) - U c fb+b n c+ n c an } + ^2 iK(b+b n 

n 

+ J2tinUn-Nn f), (10) 

a 

where all variables c+„, c an , /+„,/<*„ , b+, b n and A„ are 
functions of the imaginary time r. Let us use the following 
transformation : 



exp(U cf ^2 J b+b n c+ n c an dT) = 
an 

Ucfir- 1 J D* n D** n exp(-C/ c/ ^ J (*nK 

n Q 

+V n b+b n + C anCan)dT, (11) 

a 

where & n (T) and \&£(t) are complex conjugate Bose 
variables. Then, substituting this transformation into 
Eq. (|l0|) , we obtain the partition function as a path in- 
tegral over the Grassmann variables c+„, c an , f^ n ,fan 
and Bose variables 6+, b n , "J, \&* and A„. The integra- 
tion over the Bose variables can be performed using the 
saddle-point method that gives an exact solution of the 
model with the Hamiltonian (^J) in the limit N — ► oo. We 
look for an uniform saddle-point solution: 



«A„(t) = A, 6„(t) = b, 



(12) 



Here it should be noted that in the saddle point the \f is 
not complex conjugate to 'J, i.e. \P* ^ because Re^P 
and Im^J must be considered as independent variables. 
The physical meaning of b , \& and \& becomes clear from 
the mean field equations: 



b=- 



V 



(A + U ef 9) N u ^ 

N J ak 



( C ok/° k ) 



MF ' 



(13) 
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MF 



-Nr. 



(14) 



E f - U cf N c 



(22) 



* = - \H 



(15) 



where N u is the number of unit cells in the lattice, N c 
is the number of conduction band holes per unit cell , 
< ... >mf means averaging over the mean- field Hamil- 
tonian Hmf that easily follows from Eqs.(^0|)-(^5|): 

H M F = Yj{( £k ~ Uc S H 2 ) C ak C "k + £//« k /ak 
ak 

+(Vb*c+ k U + H.c.)} - N u (e f - Ef)Nf, (16) 



where et is the renormalized / level energy: 
e f = E f + X, 



(17) 



N f = J2 a (fanfan) is the average number of / holes on 
the / shell on the site n. The parameter determines the 
valence (v) of the / ion. For example, for an Yb ion we 
have v = 2 + Nf. Below, for the sake of simplicity we shall 
assume that there is only one Yb ion per unit cell. The 
order parameter b describes the formation of the HF state 
(in other words, the coherent Kondo state). According 
to the constraint Eq. (|2j) (or Eq.(||)) there is the following 
relation between b and a valence change (AN/) of / ions: 

\b\ 2 = Nn of - J2 (fLUn) = No/ -N/ = ANf (18) 

a 

Because in the HF state |6| ^ 0, we come to the very 
important conclusion that in the HF state / ions have a 
non-integral valence. For example, for Yb ions we have 
Nq/ = 1, and the valence is equal to v — 3— AiV/. Strictly 
speaking, due to valence fluctuations governed by the hy- 
bridization terms in the Anderson Hamiltonian Eq. ([!]), 
even in the high temperature region where (b n ) = the 
valency is not integral, i.e. N/ = Nq/ — (b+b n ) < N /, 
because {b^b n ) ^ 0. In order to take into account the ef- 
fect it is necessary to go beyond the mean-field approach 
and calculate 1/N corrections. 

It is convenient instead of the dimensionless order pa- 
rameter b to use the following energy parameter charac- 
terizing the effective hybridization: 



B = V*b 



(19) 



Inserting Eq.([L9|) into Eq.([L8|) gives a useful relation: 

\B\ 2 = \V\ 2 AN f . (20) 
Then Eq.([j"3|) takes the conventional form: 

ak 

where the exchange interaction J is given by: 



From Eq.(|l6|) it follows that the interaction U c / leads to 
the energy shift U c fANf of the conduction band to lower 
energies. This in turn shifts the chemical potential [i and 
effective energy s/ with respect to their bare values fio 
and £o/. 

In accordance to Eq. (|2^) the interaction U c f influences 
the exchange interaction J, which in turn effects on the 
Kondo screening. This effect plays a very important role 
and will be discussed below in detail. From the physical 
point of view the effect is related to the renormalization 
of the energy of the first excited state of Yb ions due to 
the Coulomb interaction. 

The diagonalization of the Hamiltonian Hmf shows 
that in the HF state the renormalized band structure 
consists of two hybrid bands: 

^k = i{£k + e f T [(ek - e/f + 4 |S| 2 ] 1/2 }- (23) 

where the upper and lower signs corresponds to 77 = 1 
and 2, respectively, and the following energy parameters 
are introduced: 



Jf=£f + U c /AN/, 

E vk = £„ k + U cf AN f , 
Ji = H + U c fANf. 



(24) 



The free energy of the system is determined by the well 
known equation: 



F = E- ST 



(25) 



The internal energy E and entropy S of the electron sys- 
tem per unit cell (or per / ion, because we assume that 
there is only one / ion per unit cell) is equal to 



N 



E = ^EE *W(^k) - U c f AN f N t (26) 

u ?)=1,2 k 

-(e f -Ef)Nf, 



iVu 77=1,2 k 

+(1-/(^0) Mi-/^))}, 



(27) 



where 



N t = Nc + Nf = N 0c + No/ 



(28) 



is the total number of holes per unit cell, and f(E) = 
(exp((E - /£)/T) + I)" 1 = f(E). It is interesting to 
note that iV c = Nq c + ANf, which demonstrates the fact 
that a certain part of holes from / shells transit into the 
conduction band. 
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III. MEAN-FIELD EQUATIONS 

The state of the system at any temperature T is com- 
pletely determined by three parameters: ju, B and £/. To 
determine the parameters it is necessary to solve a set of 
three nonlinear mean-field equations: 



iV * = ^E(( C ^k)+(/ Q + k /« k )): 
ok 



Nt 



B 



Y. ( C ak/«k) ■ 



(29) 



(30) 



(31) 



At arbitrary T the mean-field equations have a trivial 
solution B — 0. Moreover, nontrivial solutions with B ^ 
can also exists. These solutions correspond to either 
minimum or maximum of the free energy F. In the case 
B =/= one can rewrite the mean-field equations in a form: 



N t = N dep( £ )(f(E lk )+f(E 2k )), 



(32) 



N f = N dep(e) 



|5| 2 /(£ik) + (e/-£ik) 2 /(£ 2 k) 



j = nJ dep(e) 



E-2k — -Elk 



(33) 



(34) 



where p(e) is the density of states (DOS) in the conduc- 
tion band: 



(35) 



It is convenient to represent the internal energy E, 
Eq. (p6|) , as follows: 

E = NJ2[ dep(e)E vk f(E nlc ) + - \B\ 2 

V =1,2 J 

+ ^\Bf-(e f -E f )N 0f . (36) 

Let us analyze the exchange interaction J given by 
Eq.@. Using the relations N c = N 0c + AN f and 
Eq.(20), one can rewrite J 1 in the form: 

1 1 



/ 2U B 2 \ 

— = — y [Ef-Ef- U cf N 0c 2* (37) 

J \ V \ 2 I 1 ct \ V \ 2 I V ' 



Eq.(p7|) shows that in the HF state, i.e. at B ^ 0, the 
exchange interaction J is enhanced with respect to one 
in the normal state with B = 0. It is the effect that is the 
driving force of the first-order valence phase transition in 
YbInCu4 as it will be shown below. 

In the case U c f — at T < 7\. the mean-field equations 
Eqs.(^2|)-(Q) have only one nontrivial solution with B ^ 
which describes a continuous transition into the HF 
state. We call the solution "conventional". The Kondo 
temperature Tfc can be found from Eq.([34"|), supposing 
B = 0. A detailed analysis of the HF state described 
by the conventional solution can be found, for example, 
in the review. p0[ | Temperature behavior of ANf and 
susceptibility in the HF state found from a numerical 
solution of Eqs.(|3^)-(^i|) will be discussed below in Sec. 
VI. 

One can consider the formation of the Fermi-liquid 
state with a non-zero effective hybridization B as a co- 
herent Kondo effect . [^o| The mean-field theory enables 
us to find the leading contribution in terms the l/N ex- 
pansion into the phenomena. In the high temperature 
region the mean-field solution with B = discribes a 
localized / holes that do not interact with conduction 
holes. To take into account the interaction it is neces- 
sary to consider l/N corrections due to slave boson fluc- 
tuations around the mean-field solution. The corrections 
can be important in the region of crossover from a weak 
coupling regime to a strong coupling regime that occurs 
at T ~ T k . f§ 

In addition, one can note that at T = Tk the Green's 
function < b^{iuj)b n {iuo) > has a pole at uj — 0. It is the 
singularity that results in the Abrikosov-Suhl resonance. 

At low temperatures T << Tk in the HF state a new 
energy scale arises. It is the so called low temperature 
Kondo scale To determined as 

T = E f - p. (38) 

Solving Eq.(H) at T = in the case p(e) = po gives 



T = (e/ - -Emir,) exp(- 



Np J 



(39) 



where -Emin = E\(k — 0). It is important to note a dif- 
ference between the two Kondo scales T and Tk- But 
the Abrikosov-Suhl resonance, the Kondo temperature 
Tk appears also in a high-temperature expansion of the 
spin susceptibility, for example, and other physical pa- 
rameters in the weak coupling regime. To is an energy 
scale that determines thermodynamic properties in the 
Fermi-liquid regime with strong coupling. In general 
case, To ^ Tk (for a detail discussion on the problem 
see, for example, a review pq]). 



IV. STRICTION PHENOMENA 

When an Yb ion transits from the ground state Yb +3 
into the excited state Yb +2 , the radius of the / shell 



G 



increases. It results in an local internal pressure that 
brings about a local lattice strain. Let e^, where i,j = 
x, y, z, be a local strain tensor. Since the probability to 
find an Yb ion on the site n in the divalent state is given 
by the operator b^b n , we can write the energy related to 
this type of the electron-lattice interaction in the form 



H e -i = - ^dijeji{n)b+b n . 



(40) 



For a cubic lattice the energy tensor dij has a simple form 
dij = Sijd , therefore this electron-lattice interaction can 
be rewritten as: 

ffe-i = -dJ2eB(n)b+b n . (41) 

n 

where e B i n ) = g xx (n) + e yy (n) + e zz (n) is a local volume 
strain. For Yb compounds the interaction energy d is 
positive, which corresponds to a local volume expansion 
when Yb ions change their valence from +3 to +2. In 
turn, d is negative for Ce ions with the ground state 
configuration 4/ 1 , because the / shell of a Ce ion in the 
excited 4/° state has a smaller radius than in the 4/ 1 
state. In order to make clear a physical meaning of the 
interaction H e _i let us find an average value < H e _i > 
in the case of an uniform strain e b (h) = e B ■ Taking into 
account Eq.(fL8J) we find 



< H e -i >= -de B AN f N u 



(42) 



The volume strain es is related to the volume change 
(At?) of the whole volume (v), namely, eb — Av/v. Then 
Eq.([42|) takes the form 

< H e -i >= poANfAv (43) 

where we introduce a new fundamental parameter 



pa = -d/v . 



(44) 



Here vq is the unit cell volume (per / ion) . The parameter 
Po has the following physical meaning: pg is a pressure 
produced by a / ion on the lattice when the / ion valence 
is changed by 1. po is negative for Yb ions, because it 
tends to expand the lattice, and positive for Ce ions, 
because it tends to shrink the lattice. The energy Eq. ( |43| ) 
has a simple physical meaning. This is the work produced 
by the pressure poANf in order to change the volume of 
the system by value Av. 

Apart of the /-shell-size-fluctuation mechanism 
(Eq.(^Oj)) of the electron-lattice interaction there is an- 
other mechanism related to a volume dependence of the 
hybridization parameter V in Eq. (Jl|) : 



V{e B ) = Vexp(re B ) 



(45) 



where r < that corresponds to increasing Vies) when 
decreasing volume. Eq.([45|) follows from a volume de- 
pendence of overlapping between / and conduction band 



wave functions. Usually it is supposed that it is the 
mechanism that is responsible for a pressure dependence 
of Tfe and the Kondo volume collapse phenomena in 
compounds with Kondo impurities and HF compounds. 
However, as we shall show below it is not valid in the case 
of Yb compounds in which the /-shell-size-fluctuation 
mechanism plays a more important role. Taking into 
account these two mechanisms of the electron-lattice in- 
teraction represented by Eqs.(|4l|) and fliq), we obtain the 
following total Hamiltonian: 



H = H A + H C + fle-l + Hia 



where 



Hi 



1 



1 a I 



-G B e B N u 



(46) 



(47) 



is the energy of a uniformly deformed lattice. Cb is the 
bulk modulus per / ion. The Hamiltonian Eq. (^6|) results 
in the same mean- field equations (|3^)-(|34|) with the only 
difference. Namely, the exchange interaction J becomes 
strain dependent: 



J{e B ) 



\V{e B )Y 



If-Ef- U cf N 0c - 2U cf AN f - de B 



(48) 



Taking into account the elastic energy Eq.(|47j) and elec- 
tron energy Eq. (|36|) , the total internal energy per / ion 
can be written in the form 



E t = E+^C B e 2 B . 



(49) 



In the equilibrium state the strain es is determined by 
a minimization of the free energy with respect to e# : 
dF/de B = 0. Taking into account Eqs.® and @ the 
minimization results in an equation: 



-A 

■C B 



2r\V(e B )\ 



C B J(e B ) 

which easily follows from the relation 



)AN f 



(50) 



dF/de B = dE t /de B 
d 



\B 



(- 



de B J(e B 
Cb&b 



' ) + U cf \B\ i 



d 



de B \V(e B )\ 



-) 



since the entropy S, Eq . (|27j) , does not depend on es in 
a direct way (here we neg lect a volume dependence of 
the Fermi energy). Eq. (|50|) shows that the volume strain 
es is proportional to the valence change ANf of the Yb 
ions. Such a proportionality have been observed experi- 
mentally. fll| Eq.@ together with Eqs.(||)-(|3|), where 
J is replaced by J(e B ), represent a closed set of nonlinear 
equations which determine the equilibrium parameters 
e>(T), Ji(T), B(T) and e B (T). 

In the case of a sufficiently small strain |es| << 1 one 
can use a linear approximation: 
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(51) 



where po is the DOS on the Fermi surface, f2 = fix + ^2 " 

n 3 , 



fii 



9 In J(e B ) 
c9 In c 



(52) 



J3.Ef.CB —0 



2r 



2 = — 



«3 = 



1 fdJ- 1 dB dJ- 1 de f 



Np \ dB de B def des 



Since r < 0, consequently, f2i is always positive while 
1^2 can be both positive and negative. In general case f2 
can be both positive and negative. The physical meaning 
of the parameter comes from Eq. ([j9|) . Substitution of 
Eq.Q into Eq.@ gives 



T (e B ) ~ exp(- 



1 



N Po J(0) 



fie B ). 



(53) 



The equation determines the dependence of the low tem- 
perature Kondo scale To on the strain e B . 

Electron-lattice interactions described by Eqs.([h)]) and 
( [i"5| ) lead to a renormalization of the elastic constants. 
In order to find the effect it is necessary to calculate 
the free energy F as a function of e B at given T : F = 
F(T, e B ). For this purpose we have to solve Eqs.(H|)-(|3j) 
with the coupling Eq. (|4^) at fixed e b ■ This enables us to 
find £/(T, e B ), iu(T, e B ) and B(T, e B ) and then F(T, e B ). 
Renormalized elastic constants c* 7 - are equal to 



(Zeeman energy). The Zeeman energy of the spins in a 
magnetic field H directed along z axis is equal to 



H z = -gp B Hj2(S z fn + S z cn ), 



(56) 



where g and \ib are the gyromagnetic factor and the Bohr 
magneton, respectively, S* n and S z j n are z components of 
the spin operators for conduction band and / holes: 



E 

a— — 

3 

Sfn — E a fanf< 



a=— 3 
3 

" „,f+ t 

can • 



(57) 



The magnetic field splits the TV-fold degenerate hybrid 
bands Eq.(pl): 



Enak — E„ k — agpsH. 



(58) 



Neglecting a magnetic field effect on the orbital motion of 
heavy fermions, which is possible due to the large mass of 
the quasiparticles, one can write the mean-field equations 
(§§)-(§f) in the form: 



^ = E/ dep(e)(f(E lak ) + f(E 2ak )), 



(59) 



N f = E/ d£ P^ 



\B\ 2 f{E lak ) + (s f - E lk ff{E 2ak ) 



\B\' + (e~ f -E lk y 



(60) 



d 2 F 



(54) 



33 / e=e(T) 



where e^(T) = es(T)/3 is the equilibrium elastic strain. 

An effect of applied pressure P on the system under 
consideration can be calculated from the thermodynamic 
relation 



P= - 



dF{v,T) 
dv 



(55) 



Eq.(p5[) together with Eqs.(|32|)-(|3^), where J is replaced 
by J(es), represent a closed set of nonlinear equations 
that determine the equilibrium parameters £/(T, P), 
ju(T,P), B(T,P) and e B (T,P). 



V. MAGNETIC FIELD EFFECT 

Let us study an influence of a magnetic field H on 
the system and the valence phase transition. We shall 
take into account only an interaction of spins of conduc- 
tion band and localized / holes with the magnetic field 



~ = Y [ dsp(e y f{El Z k) ~ L {E2ak) . (61) 
J E 2V -E, 



If we take into account the electron-lattice interaction 
then the exchange interaction J in Eq.(|6"l|) is given by 
Eq. (ff8|) . Solving these equations together with Eq. (pT)|) , 
we can find parameters e/(T,H), p(T,H), B(T,H) and 
e B (T,H). Then we can calculate the total internal en- 
ergy per / ion: 

Et = y y I dep(e)E vak f(E vak ) + — ^— \B\ 2 

a 77=1,2 J[ - eB > 

+ , T ^ C/ „4 l^l 4 - - E f )N of + \c B e\, (62) 
\V{e B )\ 2 

and the entropy 

S = ~y y f dep(e){f(E nak ) In f(E vak ) 

a 77=1,2'' 

+ (1 - f(E vak )) ln(l - f(E vak ))}. (63) 
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Calculating F as a function of T, H and e# enables us 
to find the influence of the magnetic field on the elastic 
constants Eq.(Q). 

In the heavy-fermion state the magnetization of the 
system per unit cell can be given as either a sum over 
magnetic moments of conduction band and / holes or a 
sum over magnetic moments of heavy quasiparticles in 
the hybrid bands Eq.(|58|): 



M 



MBN- 1 J2(< S cn > + <Sf n >) 



9fJ>B 



a 7)=1,2 



dep(e)f(E v 

nk ) • 



(64) 



Differentiating M with respect to H at H = one can 
find the static magnetic susceptibility per / ion: 



X(T) 



-1(9Vb) 2 3(3 + 1)N J2 J dep(e)f(E vk ), 



77=1.2 ' 



(65) 



where f'(E) = df(E)/dE. Moreover, we took into 
account that at small magnetic fields the parameters 
e f (T, H), p(T, H), B(T, H) and e B (T, H) have field cor- 
rections of order of 0(H 2 ). At T = we obtain the well 
known result (see, for example, the review p0[|): 

X (T = 0,H = 0) = \{gp B ?j{j + l)Np* F . (66) 

Here the renormalized density of state p* F and the effec- 
tive quasiparticle mass to* on the Fermi surface are given 
by 



Pf 
Po 



UL = = 1 

m 



\B£ 



(67) 



In the right hand side of Eq.(|67|) the first term is con- 
tributed by conduction holes, and the second term is con- 
tributed by / holes. ^,^8| It is obviously that in the case 
rn* /mo >> f the zero temperature susceptibility (^) is 
mainly determined by / holes. In the case po » To at 
T = H = from Eqs.(p2[)-(p3|) we obtain a useful relation 



\B 



,2 = N f To 
N Po ' 



(68) 



An additional assumption Np T << 1 gives the follow- 
ing result: J|(| 

X (T = 0,H = 0) = \{gp B fj{j + 1)^. (69) 

However, the inequality po >> Tq can be invalid for a 
semimetal similar to YbInCu.4 having a low charge carrier 
concentration and large enough To. Therefore in general 
case it is better to use Eqs.(|6oj) and (pTj). If T >> 
T v where T v is the critical temperature of the valence 



phase transition, then at T < T v in the HF state the 
susceptibility x(T, H — 0) will have a weak temperature 
dependence, as a temperature correction to x(0) is of 
order of O(T 2 /T 2 ). 

Above T v , that is in the normal paramagnetic state 
with incoherent Kondo scattering, the total susceptibility 
of the electron system considered is mainly determined 
by localized / spins, as at temperatures T << p the 
Pauli susceptibility of conduction holes is much smaller 
than the susceptibility of the localized / holes weakly in- 
teracting with conduction holes. Within the mean-field 
approach the interaction is neglected, and / holes behave 
as free paramagnetic spins with the Curie- Weiss suscep- 
tibility: 



1 



X (T,H = 0) = -(gp B yj(j+iy- 



N, 



of 
T ' 



(70) 



This result is valid in the leading order in l/N . In order 
to find logarithmic corrections into x due to the Kondo 
screening, it is necessary to study local fluctuations of 
slave bosons bi around the paramagnetic state. Q 

One can conclude that on decreasing the temperature 
below T v the susceptibility undergoes a jump from the 
value Eq.(|7(]) to much the lower value Eq.(p6|). The result 
is in qualitative agreement with experimental data (see 
paper [fiol and recent data fl9|,|22|| ) . 



VI. RESULTS OF NUMERICAL CALCULATIONS 

In the framework of the mean-field approach to the ex- 
tended lattice Anderson model at a given temperature T 
and pressure P the free energy is completely determined 
by four physical parameters: the chemical potential p, 
the effective energy Sf of the / level, the order parame- 
ter (effective hybridization) B and the volume strain e b ■ 
In order to find the parameters we solved numerically the 
set of equations (|32j)-(|34|) and (|55|) with the exchange en- 
ergy Eq.(48). For the sake of simplicity we used a fiat 
conduction band with an energy independent DOS pq. 

We found that, depending on the / level energy Et, the 
hybridization parameter V , the conduction band DOS po, 
the initial concentration Nq c of conduction band holes 
and the interaction U c f , the set of the mean-field equa- 
and (F 



tions (|32|)-(|34J) and ( |55| ) has different solutions that give 
different scenarios of temperature behavior. 

At first let us briefly discuss different solutions in the 
simplest case when the electron-lattice interaction is ne- 
glected and P = 0. At U c f = and T > T k the mean- 
field equations has only a trivial solution with B = 0. 
At T < Tk apart from the trivial solution there is one 
nontrivial solution with B =/= 0, namely, the conven- 
tional solution describing a continuous transition into the 
HF state. The temperature dependence B(T) of the con- 
ventional solution is represented in Fig.|l| at parameters 
N 0c = 0.069, \V/p \ 2 = 19 and l/Np J = 3 . At T < T k 
the HF state with the B{T) has a lower free energy than 
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the trivial solution. Solving numerically Eqs.(32)-([34|), 
we have found that at T = and U c f = within the con- 
ventional solution the valence change AN* can achieve a 
maximum value of order 0.05. 

With increasing U c f at fixed parameters Ef, V, po 
and Nq c ~ 0.07 the behavior of ANf (T) remains qual- 
itatively the same as for the conventional solution with 
U c f = 0. There are only quantitative differences which 
are nevertheless very important. At a certain value of 
U cf (U c f/(Np \V\) ~ 3) the valence change AN f (0) 
can achieve a value of order 0.1. Such a scenario we 
shall use below for describing the physical properties of 
YbAgCu 4 . 

At large enough U cf {U cf /{Np \V\ 2 ) > 5 and Y c0 ~ 
0.07) in a certain temperature region the mean-field equa- 
tions (^)-(j^) have two nontrivial solutions correspond- 
ing to minimums of the free energy F as a function of 
the effective hybridization B (there are also solutions 
corresponding to maximums of F(T,B)). At parame- 
ters U cf /(Npo\V\ 2 ) = 5.6, N Qc = 0.069, \V/p \ 2 = 19 
and 1/NpoJo = 3 (below we shall use these parameters 
for studying properties of YbInCu4) these solutions are 
represented in Fig.[l| One solution with small B begins 
at T = Tfc , i.e. B(Tk) = 0, and ends at a point s 
at temperature T s — 0.83Tfc. We shall call the solution 
"soft". Near Tk the soft solution behaves similar to the 
conventional solution. Moreover, there is another solu- 
tion with larger B. It starts at a point h at temperature 
Th ~ 1.06Tfc, exists up to T — and has a weak temper- 
ature dependence. We call the solution "hard" . The free 
energy F(T, B) as a function of B at different temper- 
atures T is represented in Fig.^. At T > Th the F has 
only one minimum at B = 0. In the range Tk < T < Th 
the F has two minimums, one at B = and the other 
at B corresponding to the hard solution. At T < Tk 
the trivial solution gives a maximum of F. In the range 
T s < T < Tk the F has two minimums related to the soft 
and hard solutions. Below T s there is only one minimum 
given by the hard solution. For the given parameters a 
first-order phase transition takes place at T v w 0.97Tfc . 
At T > T v either soft or trivial solutions give an absolute 
minimum of F , while at T < T v the absolute minimum is 
given by the hard solution. The valence jump ANf is of 
order 0.2. This scenario will be used below for describing 
the first-order valence phase transition in YbInCu4. 

At even larger U c f the order parameter B correspond- 
ing to the hard solution can be much larger than B of 
the soft solution. The hard solution can arises well above 
Tk since in this case Th » Tk- The soft solution can ex- 
ist down to very small temperatures since in this case 
T s «Tk- Besides, a first-order valence phase transition 
can occur at a critical temperature T v » Tk with the 
valence change AJVy(0) ~ 1. Such a scenario occurs, for 
example, for the model parameters U c f/(Npo \ V\ 2 ) = 10, 

N 0c = 0.069, \V/p \ 2 = 20 and l/Np J = 4 which give 
T v /T k = 130. 

Our analysis of the numerical solutions reveals that 



there is a correlation between ANt(0), T). and T v . If 
Tk << T v , which takes place at sufficiently large U c f 
(U cf /(Npo \V\ 2 ) > 9), then AN f (0) is close to 1. If T v is 
only slightly larger than T k , then AJV/(0) ~ 0.25. In the 
frame of the first-order phase transition scenario a mini- 
mum value ANf(0) ~ 0.2 is achieved when T v is slightly 
smaller than Tk- 

There are also regions in the parameter space (£7, V, 
Po, M)c, U c f) in which even at U c f 7^ a first-order phase 
transition does not occur. If the electron-lattice interac- 
tion is included then one can expect a rich T — P — H 
phase diagram. 



A. Physical properties of YblnCiu 

Let us apply the numerical calculations discussed 
above for studying the thermodynamic properties of 
YbInCu4. For the sake of simplicity we suppose the con- 
duction band to be flat and the degeneracy N = 8. In 
accordance with the Hall measurements |llj we use the 
initial charge carrier concentration iVo c =0.069 per for- 
mula unit. The best fit to the experimental data is found 
at the following parameters: V = 0.267 eV, U c f = 0.451 
eV, iV> =1.125 eV" 1 , p -E f = 0.272 eV. The param- 
eter NpgV 2 characterizing the / level broadening due to 
the hybridization is calculated to be 0.08 eV. One can 
note that for rare earth compounds the broadening is 
typically of order 0.01-0.1 eV. H The initial chemical 
potential po is calculated to be 0.061 eV. The obtained 
value of V is in good agreement with the mixing terms 
0.27 eV between Yb 4/ and the Cu p states found in the 
electronic band structure calculations [|l3| of YbInCu4. 
According to the calculations the Cu p states give 
the main contribution into the Ti states in the vicinity 
of the Fermi surface. 

The energy of the 4/ level, p - E f = 0.272 eV, is 
consistent with the experimental results of PES , 0.3 eV, 
for the compound. |32j One can note that in other Yb 
based mixed valent compounds such as YbAl2 and YbAl3 
the 4/ level energy estimated from photoemission spectra 
pl| , p2[ also is not centered at the Fermi energy, but 0.24 
eV below it. 

For describing electron-lattice interactions it is neces- 
sary to set the parameters r and d. In many HF com- 
pounds the parameter r is of order — 8 < r < — 2. We take 
r = —2 and d = 2 eV that give fii « 12 and SI2 = —25 
in Eq. ( |52| ) . The unit cell of YbInCu4 contains four Yb 
ions and has volume 372x 10 _24 cm 3 . Therefore each Yb 
ion occupies volume i>o =93 xl0~ 24 cm 3 . Since according 
to experimental data |l3|,fl4} in the normal state the bulk 
modulus of YbInCu4 is equal to 11,1 x 10 11 erg/cm 3 , 
the parameter Cb in Eq.(47) is taken to be equal to 
1,03x10" 
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c-bVq = l,Udx ill " erg. 

Our analysis of the numerical solutions of the mean- 
field equations (^)-(^) and (55) shows that for the pa- 
rameters chosen the system under consideration under- 
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goes a first-order isostructural valence phase transition at 
the critical temperature T v =42 K. Above T v the system 
is in the normal state. Below T v the HF state described 
by the hard solution is formed, as the free energy of the 
HF state becomes smaller than the free energy of the 
normal state with the incoherent Kondo scattering. It 
is demonstrated by Fig.|[ Within the mean-field solu- 
tion at T > T v the Yb ions have the integral valency +3 
corresponding to B = 0. At T < T v the Yb ions have 
a non-integral valence 3 — ANf . A temperature depen- 
dence of the ANf is shown in Fig.^. At T = T v there is a 
jump of the valence, ANf = 0.24. With decreasing tem- 
perature the valence change ANf achieves the value 0.27. 
The value is larger than 0.17 deduced from L3 measure- 
ments. jll| Reasons of the discrepancy will be discussed 
below. 



1. Thermal properties of YbInCv,4 

In FigjH|we represent the temperature behavior of the 
entropy (5) of the system considered. In the framework 
of our approach above T v the entropy is large enough 
due to a large contribution given by localized / holes. At 
T v the entropy drops from S(T V + 0) = 25 J/mol K to 
S(T V - 0) = 4.9 J/mol K, i.e. AS = 20 J/mol K. These 
theoretical estimates are larger than experimental data: 
AS = 10 J/mol K,S(T v +0) = 13 J/mol K, S(T v -0) = 3 
J/mol K. |l{| However, the theoretical estimate of the 
ratio S (T v + 0) / S (T v — 0) = 5.1 is in good agreement with 
the experimental value S{T V + 0) /S(T V - 0) = 4.3. Below 
T v f holes are strongly hybridized with conduction band 
holes and form the nonmagnetic heavy fermion Fermi- 
liquid ground state. That is why the entropy tends to 
zero in the limit T — » 0. 

There are a few reasons of the divergences between 
our estimates and experimental data.. At first, we ne- 
glected a crystal field splitting of the 8-fold degenerate 
/ level. Inelastic neutron measurements ji5| revealed 
that the quartet Tg, is the ground state, but the split- 
ting A c is quite small, A c = 32 K. The experimental 
value S(T V + 0) = 13 J/mol K« i?ln5 is close to the 
value i?ln4 expected for four-fold degenerate Tg state. 
An account of the crystal-field splitting within our ap- 
proach will result in decreasing entropy and a more re- 
markable temperature dependence of the entropy above 
T v . At second, the mean field approach does not give a 
correct result of the entropy of the localized / electrons 
in the high temperature phase. Namely, the approach 
results in the value S = R(lnN — (N — 1) ln(l - A^ 1 )) 
for the entropy of noninteracting / holes instead of the 
value i?ln N. To obtain a correct result it is necessary to 
go beyond the mean-field solution and take into account 
fluctuation corrections in next orders in 1/N. 

In the HF phase at T < < T v the linear coefficient of the 
specific heat 7 is given by the relation 7 = hir 2 k%N p* F 
where p* F is determined by Eq. (|67j) . For the parameters 



chosen we find 7 = 71 mJ/mol K 2 . The experiment |Tg| ] 
gives slightly smaller value 7 w 50 m J/mol K 2 . 

2. Magnetic susceptibility of YblnCui 

Using Eq.(|65|) we calculated the temperature depen- 
dence of the spin susceptibility x(T). Corresponding 
results are represented in Fig.||. Above T v localized / 
electrons give the main contribution into the suscepti- 
bility. As a result, the x follows the Curie- Weiss law: 
X = C/T. Below T v hybrid quasiparticles are formed, 
and the electron system under consideration behaves as 
a Fermi-liquid system with an enhanced Pauli suscepti- 
bility. At zero temperature the x(0) is given by Eq.(p6|). 
At ambient pressure we found T L = C/x(0) = 384 K. 
The experimental value is 470 K. jnj We calculated also 
the low temperature Kondo scale To defined by Eq.(|3^) 
and found T = 315 K. In the conventional HF sys- 
tems it is the To that plays the role of a universal en- 
ergy scale for low temperature thermodynamic and trans- 
port phenomena in the HF state (see Refs. |30| and 
|Q, for example). According to Eq.(|69[), in the limit 
fiQ » Tq there is a simple relation between To and Tl : 
T = Nf T L = (1 - AN f ) T L . The approximate rela- 
tion at ANf = 0.27 gives Tj= 430 K that is slightly 
larger than our calculated value (T L = 384 K). The dif- 
ference is a consequence of the approximate character of 
the relation, because in our analysis we have Ho/Tq ~ 2.2 
due to the small hole concentration N 0c — 0.069 per / 
ion. Our calculations revealed that both the Tl — 384 K 
and To = 315 K are much larger than the Kondo temper- 
ature Tfc = 40.4 K characterizing Kondo effect at T > T v . 
The result is in agreement to inelastic neutron scattering 
measurements |45|Ji"7| , |48[ that also revealed an enhance- 
ment of the energy scale from Tk = 25 K to To = 405 
K. The very large enhancement of To in comparison with 
Tfc is brought about by the enhancement of the exchange 
energy J. In papers Jll|,|l2| in order to explain the en- 
hancement of Tl it was suggested that the effect is caused 
by a strong energy dependence of the conduction band 
DOS p{e). Namely, p(ef) must be sharply increasing 
function of the Fermi energy e f ■ The results of our cal- 
culations show that in terms of the extended lattice An- 
derson model the enhancement can be explained even in 
the case of a flat conduction band. 

The enhancement of 7 and x(0) with respect to the 
values typical for normal metals is explained by the in- 
crease of the DOS on the Fermi surface. According to 
Eq. (|67|) , there is also a quasiparticle mass enhancement. 
Our calculations give m* /mo = 27. Therefore, one can 
classify YbInCu4 as a "light" heavy-fcrmion compound 
as it was suggested by Felner et al. 0] 
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3. Volume expansion and bulk modulus in YblnCu^ 



5. H — T phase diagram. 



In FigJ^ we represent results of the numerical calcula- 
tions of the volume strain es = Av/v. In good agreement 
with experimental data |5 11 44|] we find a sharp volume 



expansion about 0.5% when decreasing T below T v . Cal- 
culating the volume strain at different pressures, we es- 
timated the isothermal compressibility fcx = —dlnv/dP 
and found k T (T = 20 K) = 1.05 Mbar" 1 . Experimental 
measurements Hgave k T (T = 20 K) = 1.2 Mbar" 1 . Ul- 
trasonic studies |[|||i9| resulted in k T (T = 20 K) = 0.95 
Mbar" 1 . 

Using Eq.([54|) we calculated the temperature depen- 
dence of the bulk modulus cb — 1/ fcr- These results are 
presented in Fig.|[ Such a temperature behavior is in 
complete agreement with the experimental data. fl43p9| ] 

According to Eq. (|50|) , at a sufficiently small electron- 
lattice coupling the volume strain must be proportional 
to ANf. es — aANf. Our calculations conform the re- 
lation and give a=0.019. The experimental estimation of 
the linear coefficient a is 0.046. 11 



4- Pressure effect in YblnCui 

The effect of pressure on the physical properties of the 
system under consideration is represented in Fig|6| Ac- 
cording to our calculations, applying pressure shifts the 
critical temperature T v of the first- order phase transi- 
tion to lower temperatures at the rate dT v /dP = —0.13 
K kbar -1 . The theoretical estimate is one order of mag- 
nitude smaller than the experimental result ]To| , ^9| , ^4[ 
, —2.2 K kbar" 1 . In order to improve the result, it 
is necessary probably to go beyond the mean-field the- 
ory. Pressure results in decreasing the valence change 
at the rate dAN f {3.5 K)/dP = -6.3 Mbar" 1 in a sat- 
isfactory agreement with the experimental value —4.5 
Mbar" 1 found from volume expansion measurements. 
p4| In Fig.p we plot magnetic susceptibility at ambi- 
ent pressure and P = 3 kbar. From the results we find 
dT^/dP = —12.7 K/kbar in good agreement with the 
experimental value —12.4 K/kbar obtained by Sarrao et 
al. The theoretical estimate of the Gruneisen param- 
eter r = — d\nTi,/d\iiv = —31 is in excellent agreement 
with the experimental value [|l9| T = —30.6. The phys- 
ical origin of the pressure effect is related mainly to the 
influence of pressure on fluctuations of the / shell size. 
Indeed, an applied pressure decreases volume, €b < 0. 
In turn, this results in decreasing the exchange interac- 
tion J (see Eq.(|48|)) and, consequently, decreasing the 
Kondo scales Tk and To. We believe that the /-shell- 
size-fluctuation mechanism explains also the negative T 
observed in other Yb based HF compounds and their 
pressure dependent properties (see, for example, Ref. ]52] ] 
and reference therein). 



Basing on the results of Sec. V we investigated the 
effect of a magnetic field on the first-order valence phase 
transition. Only the Zeeman energy was taken into ac- 
count in the numerical calculations. In our calculations 
we took the gyromagnetic factor g = 8/7 expected for 
j = 7/2. For the sake of simplicity we neglected the 
electron-lattice interaction, as it gives only about 10% 
correction in accordance to our estimates. Solving self- 
consistently Eqs.©-© we found that magnetic field 
pushes the first-order valence phase transition to lower 
temperatures. The obtained T — H phase diagram is 
represented in Fig.^J. It is interesting to note that the 
elliptic equation 



H V (T) 
H V (T = Q) 



Ty{H) 

T V (H = 0) 



(71) 



is a good fit to the critical line of the first-order phase 
transition H V (T) given by our numerical calculations. 
Recently the relationship Eq.(f7l|) between H V (T) and 
T V (H) was inferred from magnetoresistance measure- 
ments. Q Our estimate of the critical field H V (T — 
0) = 47 T is also in a satisfactory agreement with ex- 
perimental value H V (T = 0) = 34 T. The critical line 
H V (T) divides H — T plane into the low temperature- 
low magnetic field region in which the system is in the 
mixed-valence heavy fermion ground state, and the high 
temperature-high magnetic field region in which the sys- 
tem is in a normal state with stable magnetic moments 
of Yb ions. 



6. Electrical resistivity and the Hall effect in YblnCu^ 

In accordance to the measurements |l4|,^2| the elec- 
trical resistivity R(T) in YbInCu4 demonstrate two pe- 
culiarities: (i) a large drop at T = T v and (ii) a weak 
temperature dependence at T < T v . Analyzing the data 
at T > T v , one can distinguish two main contributions 
into the resistivity: a large residual resistivity Rq and a 
linear temperature dependent contribution given by elec- 
tron scattering off phonons. At T = T v the phonon 
contribution is much smaller than Rq and as a result 
R(T V + 0) ss Rq. In the low temperature phase the resid- 
ual resistivity Rq is approximately 12 times smaller than 

Rq. 

Basing on our approach let us analyze the tempera- 
ture behavior of the resistivity at T < T v . At first, one 
can note that in the HF state the renormalized electron- 
phonon interaction is very small and its contribution 
into resistivity can be neglected (see, for example, Ref. 
J39|). Only impurity scattering and collisions between 
heavy fermions determine the resistivity. It leads to 



A*T , where the coefficient A* 



5 46 1 . Since Tq » T v , the temperature dependence 
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of R*(T) must be weak eno ugh, w hich is in agreement 
with the experimental data. [[l4[|22| l 

Let us estimate and compare the residual resistivity 
above and below T v . At T > T v the residual conductivity 
can be estimated using the well known relation 



00 



e 2 poFvl F T 



e N 0c t 



(72) 



where vq F = pof/^o is the Fermi velocity, pqf is the 
Fermi momentum, pqf ~ Pof^o is the conduction band 
DOS, Tq" 1 is a rate of a potential scattering of charge 
carriers off impurities and lattice imperfections. Here we 
also used that Nq c ~ p^ F . In the HF state the resid- 
ual conductivity is determined by a potential scatter- 
ing of heavy fermions and is given by a similar equa- 
tion with the replacement of pq Fi vqf, Pof and tq by 
the renormalized parameters p* F ~ pftttl* , vf = Pf/tu* 
and r* = T m*/m . |4^j5^] Moreover, it is necessary 
to take into account that the renormalized Fermi sur- 
face determined by the total number of charge carriers 
N t = N 0c + N of , that is {{pf/pof) 3 = N/N 0c . Then, we 
find the residual conductivity at T < T v : 



1*1 \i 
e p F (v F ) 



2 N t T . 



(73) 



Therefore, we have Rq/Rq ~ N t /No c . Above we have 
found N t = 1.07 while Y Oc =0.07. It gives the crude 
theoretical estimation Rq/Rq ~ 15 in agreement with 
the experimental value p2], Rq/Rq ~ 12. Thus, we can 
conclude that the drop of the resistivity and change of 
its temperature behavior when transiting into the low 
temperature mixed-valent state in YbInCu4 are mainly 
related to (i) the decrease of the residual resistivity due 
to increase of the charge carrier concentration and (ii) 
the suppression of the electron-phonon scattering. 

Basing on our approach one can explain also the jump 
of the charge carrier concentration at T = T v observed 
in the Hall measurements. |ll],[l^] Indeed, at T > T v f 
electrons are localized and only conduction band holes 
with the concentration Nq c participate in transport phe- 
nomena and the Hall effect. At T < T v the / electron 
states are hybridized with conduction band states, which 
results in the formation of hybrid quasiparticles (heavy 
fermions). The total number of the quasiparticles under 
the renormalized Fermi surface is equal to Nt > Nq c . 
All these quasiparticles give a contribution into the Hall 
constant. 



B. Physical properties of YbAgCu4 

As we have discussed above, with decreasing the in- 
teraction U c f the system under consideration reveals a 
crossover from the first-order valence phase transition 
to a continuous formation of the HF state. We use 
such a scenario for describing thermodynamic proper- 
ties of YbAgCu4. In many respects YbAgCu4 is close 
to YbInCu4, but YbAgCu4 has a normal metallic charge 



carrier concentration. For our calculations we take 
Aoc=0.21 per Yb, which corresponds to 0.84 charge car- 
rier per unit cell. The best fitting to the experimental 
data is found at the following parameters: V = 0.24 eV, 
1/Npo=0.9 eV, pq — Ef=0.25 eV which are close to ones 
for YbInCu4. Only the interaction U c f = 0.276 eV is 
taken almost twice smaller than U c f in YbInCu4. The 
small value of U c t can be explained by a stronger charge 
screening due to a larger charge carrier concentration in 
comparison to YbInCu4. For these parameters we have 
N pqV 2 =0.064 eV, = 0.19 eV. The value of 4/ level 
energy is consistent with the result of PES |p0| , pl[ , 0.3 
eV. 

For describing the electron-lattice interactions in 
YbAgCu4 we use the same r and d as for YbInCu4, i.e. 
r = -2, d = 2 eV, and C B = 1,03 x 10" 10 erg. For 
the parameters chosen the equations (j32"|)-(|34|) and (|55| ) 
have only one nontrivial (conventional-like) solution be- 
low Tfc = 86 K. Results of numerical calculations based 
on the solution are represented in Figs.&@. 

Let us compare temperature behavior of YbAgCu4 and 
YbInCu4. In accordance to our numerical calculations 
the HF state in YbAgCu4 is continuously formed below 
Xfc = 86 K. This temperature correlates with the tem- 
perature of the resistance maximum, T max « 75 — 90 K. 
P|| According to our calculations represented in FigJJ, 
at T — ► the valence change ANf tends to 0.13 in very 
good agreement with the L3 data |Tl[] . 

Temperature behavior of the entropy is represented 
in Fig.^j. Unfortunately, entropy measurements in 
YbAgCu4 are unknown for us to be compared with our 
calculations. The calculated value of the linear coefficient 
of the specific heat 7 = 230 mJ/mol K 2 is between the 
data (UJ 7 = 220 mJ/mol K 2 , and 7 = 250 mJ/mol K 2 
found from other measurements. p2| 

The calculated magnetic susceptibility x(T) repre- 
sented in Fig.|| has a broad maximum at T m 40 K and 
its maximum value is approximately twice smaller than 
x(T — T v ) in YbInCu4, which is in good agreement with 
the magnetic measurements. [ p2[ The calculated value 
Tl = 120 K is slightly smaller than the experimental 
data (nJH T L = 150 K. 

The calculated mass enhancement in YbAgCu4 is equal 
to m* /mo = 87 . This value is in three times larger then 
the mass enhancement in YbInCu4. Therefore, one can 
classify YbAgCu4 as a "moderate" heavy-fermion com- 
pound. 

The temperature dependence of the volume strain e b 
is shown in Fig.|?]. In YbAgCu4 the zero temperature 
value es(T = 0) = 0.24 % is approximately twice smaller 
than the es(0) in YbInCu4. We find that the relation 
e B — aANf holds also in YbAgCu 4 with a = 0.018. 

Using the parameters given above we calculated the 
temperature behavior of the bulk modulus for two con- 
centrations of charge carriers, Ao c =0.21 and No c =0.Q7. 
In the both cases the calculations revealed a minimum 
of the bulk modulus. These results are represented in 
Fig. 8. One can see that with increasing the charge car- 
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rier concentration the minimum becomes smaller and less 
expressive. The results of our calculation for Ao c =0.07 
can be related to the behavior of the bulk modulus of 
YbIni_ :c Ag :r Cu4 at x « 0.3, (49] because in the case the 
charge carrier concentration is still small but the system 
already shows a continuous transition into the HF state. 
p2| The calculated behavior of the bulk modulus rep- 
resented in FigJ| is in good agreement, both qualitative 
and quantitative, with the acoustic measurements. E{t] 



VII. DISCUSSION AND CONCLUSIONS 

In the present section we shall discuss the mean- 
field approximation used in our paper and other prob- 
lems related to physical phenomena in YbIni_ x Aga;Cu4, 
and summarize results of our numerical calculations per- 
formed in the framework of the extended lattice Ander- 
son model proposed in the present paper. 

Our consideration of the model was based on the mean- 
field approximation that corresponds to taking into ac- 
count the leading order of the 1 /N expansion. The mean- 
field solution is the exact solution of the extended lattice 
Anderson model in the limit N — > oo. In the case under 
consideration the degeneracy N of the electron state 4/ 13 
is finite and equal to 8. Therefore, we face the problem of 
fluctuations around the mean-field solution. The fluctua- 
tions give 1 /N corrections to the solution. Unfortunately, 
so far there is no a detail analysis of the corrections into 
thermodynamic properties of the lattice Anderson model. 
There is only an exact solution of the single-impurity 
Anderson model with Hamiltonian Eq.([l]) for arbitrary 
degeneracy N. Q In the paper |}7j it has been shown 
that in the Fermi- liquid regime at T « To the magnetic 
susceptibility and heat capacity given by the mean-field 
solution are in a remarkably good quantitative agreement 
with the exact solution even for N = 2. For N = 8 the 
region of a quantitative agreement becomes broader. A 
noticeable error occurs only around the crossover from 
weak to strong coupling (T ~ Tk) where the critical fluc- 
tuations of the slave-boson field about the mean-field so- 
lution are greatest. However, a qualitative agreement 
takes place even in the region. With increasing temper- 
ature above Tk the mean-field solution converges quickly 
to the exact solution, as both solutions describe a para- 
magnetic state with localized 4/ electrons that weakly 
interact with conduction band electrons. As shown in 
p7| at high temperatures, the impurity susceptibility ap- 
proaches Curie behavior with logarithmic corrections of 
order 0(1/N) that are produced by slave-boson fluctua- 
tions around the normal state. 

There is the only significant defect of the mean-field 
solution related to the fact that in the limit TV — > oo the 
crossover from weak to strong coupling sharpens into a 
second order phase transition at Tk. Besides, at T = Tk 
the susceptibility has a break in the temperature depen- 
dence. 



Basing on the results let us consider the extended lat- 
tice Anderson model. In the framework of the mean-field 
approach both the single-impurity and lattice Anderson 
models are described by the same order parameter (see 
Eq.(|l9|)). However, unlike the single impurity model, the 
lattice model contains spatial correlations of the order 
parameter. It is well known from the theory of critical 
phenomena that with increasing the dimensionality, or- 
der parameter fluctuations become weaker. Therefore, 
one can expect that critical fluctuations in the 3D lattice 
model around the mean-field solution are not stronger 
than the fluctuations in the single-impurity model stud- 
ied in the paper | p7| . 

From this point of view let us discuss our mean-field so- 
lutions for compounds YbInCu4 and YbAgCu4. Accord- 
ing to the experimental data [ pd|j4^ , ^ , ^q| and our cal- 
culations, the low temperature Kondo scale To is about 
400K in the low temperature phase of YbInCu4, that 
is below T v = 42 K. As T v « To, one can expect 
that at T < T v « T fluctuations around the mean- 
field solution are weak and give small corrections with 
respect to small parameters 1/N and T/Tq. We think 
that it explains a good agreement between experimen- 
tal data and our numerical calculations of susceptibil- 
ity, entropy, valence change of Yb ions, volume change, 
bulk modulus and other physical parameters presented 
in Sec. VIA. At temperatures T > T v we are in the 
regime T > T v > Tk = 25 K, therefore, we again expect 
that corrections of order 0(1 /N) due to slave-boson fluc- 
tuations around the normal state are small and become 
noticeable only at T close to T v . As at these tempera- 
tures there are no correlations between Kondo scattering 
of electrons off different localized / electrons, 1/N cor- 
rections to thermodynamic properties are the same as in 
the single-impurity Anderson model. 

In YbAgCu4 we have a similar situation. We expect 
that at low temperatures T << To — 150 K the mean- 
field solution gives a good description of thermodynamic 
properties. The maximum error occurs in the crossover 
region with Tk = 86 K. It is interesting to note that 
according to our numerical calculations, in the lattice 
case the break of temperature behavior of x(T) at Tk, 
which is an artifact of the mean-field approximation, is 
very small in comparison to the mean-field solution of 
the single-impurity Anderson model f37f . 

In order to check the mean-field solution of the lat- 
tice Anderson model Eq. ([!]) one could use the dynamical 
mean-field theory that at the present time attracts much 
attention. [B8| For the purpose in terms of the theory it 
is necessary to solve the model Eq.(0) at different degen- 
eracy N = 2,4, 6, 8... , and then to compare the solutions 
with the mean-field solution based on the 1/N expan- 
sion as it has been done for the single-impurity Anderson 
model |37|. Unfortunately this problem is still open and 
is out of the scope of the present paper. 

Above we have mentioned that in accordance 
to the Hall measurements in compounds 

YbIni_ 2; Ag a ;Cu4 at x < 0.2 the charge carrier concen- 
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tration is low, 0.07 per formula unit. However, neutron 
scattering [^7j and susceptibility measurements Q can 
be well interpreted in terms of the single-impurity An- 
derson model with a metallic electron concentration. It 
enables us to conclude that in the compounds the charge 
carrier concentration is larger than that concentration at 
which it would be necessary to take into account Nozieres 
exhaustion principle. In the case x > 0.3 the charge 
carrier concentration quickly achieves a normal metallic 
concentration, and the problem of a low charge carrier 
concentration already does not confront before us. 

The L3 measurements |Tl| ] in YbInCu4 give a clear evi- 
dence for a non-integral valence of Yb ions in the normal 
state. However there is a profound difference between 
the origin of the non-integral valence of Yb ions above 
and below T v . In the high temperature region (T > T/.) 
the non-integral valence of Yb ions is produced by uncor- 
related electron transitions between the / shell of a Yb 
ion and conduction band. The transitions do not break- 
down the localized character of / electron states. Cer- 
tain correlations due to the Kondo effect arise only when 
the temperature lowers to temperatures about . As we 
have shown above within our approach, at T < T v the 
mixed valence state of Yb ions is caused by the effective 
hybridization between 4/ and conduction band states, 
which results in a renormalization of the conduction band 
and the formation of hybrid quasiparticle states due to 
the coherent Kondo effect. In other words, one can say 
that at T < T v the electron transitions become strongly 
correlated. From this point of view the change of the 
Yb ion valence is related to the change of the spatial 
electron distribution over electron states of all ions, in- 
cluding Cu and In ions, participating in the formation 
of the renormalized conduction band, which is consistent 
with the NQR and Knight shift measurements [|l5] p"8| |. 
This physical picture is also consistent with the theory 
of the mixed-valent state formation for a single Kondo 
impurity in terms of the single impurity Anderson model 
(see, for example, [ ^Opl]] ). 

There is another interesting problem related to the 
crystal field splitting of the degenerate 4/ level in 
YbInCii4. The inelastic neutron scattering measure- 
ments |^5| have revealed a crystal-field (CF) splitting of 
order 32 K in high temperature phase. Below T v the 
CF excitations disappear. Recently the effect of the hy- 
bridization on the CF splitting was also observed in cer- 
tain Ce based compounds of type ReNi by use of neu- 
tron spectroscopy |59|] at temperatures below Tk in the 
Fermi-liquid regime. Therefore, one can suppose that the 
hybridization and formation of the heavy fermion ground 
state at T < T v also influence the temperature behavior 
of the CF excitations in YbInCu4. 

The theory of the first order valence phase transition 
developed in the present paper can be applied also for 
studying valence transitions in Ce based and other rare 
earth compounds. Specifically, recent measurements ]6C| ] 
on Ce alloyed with 7 at.% Sc to stabilize the 7— Ce phase 
against (3— Ce formation demonstrated a similar physi- 



cal behavior at the 7 — a isostructural first-order phase 
transition as the behavior observed in YbInCu4, namely, 
localized / electrons in 7— Ce and delocalized ones in 
a— Ce, a valence change ANf w 0.2, a large enhance- 
ment of the Kondo scale from Tk ~ 60 K in 7— Ce to 
To ~ 1800 K in a— Ce and so on. However, according 
to experimental data, the electron-lattice coupling in Ce 
based compounds is stronger than this coupling in Yb 
compounds and can be a driving force of the transition. 

110 

In summary, we have used the extended lattice Ander- 
son model for investigating the first-order valence phase 
transition in YbInCu4. The model takes into account 
Coulomb repulsion U c f between / and conduction band 
holes, and two mechanisms of the electron-lattice interac- 
tion. For the model we have developed a mean-field ap- 
proach based on the 1/N expansion method. Within the 
mean-field approach we have found that the system un- 
der consideration undergoes a first-order phase transition 
from the normal state into the heavy-fermion state with 
mixed valent / ions. We have found that the Coulomb 
interaction U c f strongly influences on the exchange inter- 
action J between spins of / and conduction band elec- 
trons. The U c f enhances J and, as a result, the coherent 
Kondo effect is also enhanced. It is the driving force 
of the first-order valence phase transition. Basing on the 
model, we have carried out numerical calculations of tem- 
perature and pressure dependences of thermal, magnetic 
and elastic properties of YbInCu4 and obtained a good 
agreement with experimental results. We have studied 
the role of the electron-lattice coupling in the compound 
and found that the / shell-size-fluctuation mechanism 
of the electron-lattice interaction gives the main contri- 
bution into the volume change of the lattice and deter- 
mines the negative sign of the Gruneisen parameter. Our 
analysis have shown that the evolution of the first-order 
phase transition into a continuous transition in the se- 
ries of compounds Yblni-^Ag^Ciu can be explained by 
decreasing the interaction U c f with increasing Ag con- 
centration x due to increase of the charge carrier con- 
centration and an enhancement of charge screening. In 
the framework of the extended Anderson model we have 
found that a magnetic field pushes the valence transition 
to lower energies, and in the H — T plane the critical line 
of the first-order valence phase transition is described by 
an elliptic equation in complete agreement with experi- 
mental data 
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FIG. 1. Order parameter (the effective hybridization) B 
versus T. The conventional solution of the mean-field equa- 
tions Eqs.(||)-@ for the interaction U c f = is given by 
the dotted curve. Curves s and h represent "soft" and "hard" 
solutions of the equations, respectively, at a moderately large 
U cf (U cf = 0.451 eV, N = 8, V = 0.267 eV, Np =1.125 
eV" 1 , A r 0c =0.069 per formula unit, fi -E f = 0.272 eV) when 
Tk/no = 0.057 and T v /fio = 0.055. Here the electron-lattice 
interaction is neglected. 

FIG. 2. Free energy F of the extended lattice Anderson 
model as a function of the order parameter (the effective hy- 
bridization) B at different temperatures T/flo =0.04, 0.05, 
0.06, 0.07 for the same parameters V, N, Npo and No c as those 
in Fig.l (Tk/no = 0.057, T v /fi = 0.055). A deep minimum 
arising at temperature T/na <0.06 and B near 2 corresponds 
to the "hard" solution. The arrow shows the position of a 
minimum of the F corresponding to the "soft" solution. 

FIG. 3. Temperature dependence of the free energy F in 
the normal state (curve 1) and the heavy-fermion (HF) state 
(curve 2) in YbInCu 4 and YbAgCu 4 . In YbInCu 4 the HF state 
is described by the hard solution of the mean-field equations 
while in YbAgCu 4 it is described by the coventional solution. 
The free energies are given with the accuracy up to a certain 
energy shift. 

FIG. 4. Valence change (ANf) versus temperature T in 
YbInCu 4 (curve 1). Curve 2 shows the ANf in the case when 
the electron lattice interaction is neglected. Curve 3 is the 
ANf in YbAgCu 4 . Curve 4 is the ANf in the case when 
both the electron-lattice coupling and the on-site Coulomb 
repulsion U c f are neglected. 



FIG. 5. Temperature dependence of the entropy S(T) in 
YbInCu 4 (solid curve) and in YbAgCu 4 (dotted curve). 

FIG. 6. Magnetic susceptibility \ versus temperature T in 
YbInCu 4 at pressure P = (curve 1) and P = 3 kbar (curve 
2). Curve 3 is the \ in YbAgCu 4 at P = and T k = 86K. 
Curve 4 is the \ for the same parameters as for YbInCu 4 
but when both the electron-lattice coupling and the Coulomb 
interactions U c f are neglected. 

FIG. 7. Volume strain eg against T in YbInCu 4 (solid 
curve) and in YbAgCu 4 (dotted curve). 

FIG. 8. Bulk modulus cb versus T in YbInCu 4 (solid 
curve), YbAgCu 4 (dotted curve) and YbIno.7Ago.3Cu 4 
(dash-dotted curve). 

FIG. 9. Normalized 

critical magnetic field H V (T) / H V (T = 0) versus the normal- 
ized temperature T/T V (H — 0) in YbInCu 4 . Results of the 
numerical calculations are presented by solid squares. The 
critical line of the first-order valence phase transitions divides 
the H — T phase diagram of into two regions corresponding 
to the mixed-valence heavy fermion (HF) state and a normal 
paramagnetic (P) state with stable magnetic moments. The 
elliptic equation ( |ri| ) (dotted curve) is a good fit to the critical 
line. 
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